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Abstract 

We present a theory for the formation of a mesa shaped phonon 
pulse in superfluid ^He. Starting from the hydrodynamic equations 
of superfluid helium, we obtain the system of equations which de- 
scribe the evolution of strongly anisotropic phonon systems. Such 
systems can be created experimentally. The solution of the equations 
are simple waves, which correspond to second sound in the moving 
phonon pulse. Using these exact solutions, we describe the expansion 
of phonon pulses in superfluid helium at zero temperature. This the- 
ory gives an explanation for the mesa shape observed in the measured 
phonon angular distributions. Almost all dependencies of the mesa 
shape on the system parameters can be qualitatively understood. 

1 Introduction 

It is at present, impossible to create in superfluid helium, a constant relative 
velocity, w, between the normal and superfluid which is close to the Landau 
critical velocity. However phonon pulses can be created, which are pulses of 
normal fluid moving through the stationary superfluid with a high relative 
velocity, without superfluidity breaking down [H El [3l HI [5] . These pulses are 
tp = 10"^ - 10-^s, long. 

In experiments pQ-jS] a heater is immersed in superfluid helium at ~50 
mK. A current pulse, applied to the heater, creates a moving phonon sys- 
tem in the stationary superfluid. The initial transverse and longitudinal 
dimensions of such a phonon system is determined by the heater size and the 
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duration of the current pulse respectively. Then the phonon system starts 
to expand as it moves in the superfluid helium. There are two very different 
regimes of such evolution. At pressure above 19 bar phonons propagate al- 
most without scattering. So there is ballistic propagation. In this case, the 
phonon energy has an angular distribution with a characteristic cosine de- 
pendence. At lower pressures, the phonons strongly interact by three phonon 
processes. These are prohibited by the phonon energy-momentum conserva- 
tion laws at pressure above 19 bar. Therefore at low pressures there can be 
a hydrodynamic regime, when phonon system expands in "phonon vacuum" 
in a way similar to the ordinary expansion of a gas into a vacuum. 

In a recent paper [5] detailed measurements of the angular distributions of 
the energy in the phonon pulse were reported at zero pressure. The measure- 
ments were made at different distances from the heater, with various heater 
sizes and heater powers. The results showed a notable feature; the angular 
dependence of the energy flux has a fiat top, and concave steep sides. In 
[5] this kind of distribution was called a mesa shape. Such a phonon energy 
distributions cannot occur in the ballistic regime. In this case we would have 
cosine-like convex angular distribution, as it follows both from experiments 
at high enough pressure and from theoretical considerations (Lambert's law). 
At zero pressure we have, despite the assumptions made in the theoretical 
model in [5], a hydrodynamic regime due to three phonon processes in the 
wings of the angular distribution, where the phonon density is sufficiently 
high, as well as in the centre of the pulse. The dependencies of the mesa 
width and height on various parameters, obtained in [5], are nontrivial. So 
the hydrodynamic propagation of phonon systems in superfiuid helium is an 
important problem. 

To choose a method of theoretically analysing the phonon pulse expansion 
in superfiuid helium, we should take into account that in experiments [21 El IH 
[5] the phonon pulse duration is much greater than the time to attain a local 
equilibrium. Therefore there is a dynamic equilibrium which arises mainly 
due to scattering by three phonon processes, which have a characteristic 
time t^pp ~ 10^®s [6]. This applies to the experimental conditions [21 [3|, 
m [S]. Fast relaxation has been observed directly in the experiments with 
colliding phonon pulses [7]. There is further evidence from the fact that a 
phonon pulse propagates in superfiuid helium as a whole, with a velocity 
which is experimentally indistinguishable from the Landau critical velocity 
for phonons. For phonons with a linear energy-momentum relation, e{p) = 
cp, the Landau critical velocity is equal to c [S]. 

Thus the strong three-phonon scattering, within a pulse comprised of low 
energy phonons in liquid helium at bar, leads to a quasi equilibrium. The 
equilibrium state of the phonons can be defined in terms of a temperature 



2 



and a drift velocity, which means we can use a hydrodynamic approach to 
describe the dynamics of the phonon pulse. 

Another feature of such phonon systems is that the energy density in 
the phonon pulse is much larger than the ambient phonon energy density in 
superfluid helium. Therefore the evolution of a strongly anisotropic phonon 
system is essentially nonlinear. Nonlinear waves in superfluid helium, when 
w is small, have been studied for many years, but the analysis of nonlinear 
waves at arbitrary w has not been done until now. The first such analysis was 
made in Ref. |9j, where the phonon pulse evolution were studied by using the 
Bose-cone model for the phonon distribution function. There, the nonlinear 
equations of gas-dynamic type for the parameters of the Bose-cone model 
were obtained. These equations were solved for one dimensional transverse 
and longitudinal evolution of a phonon pulse. In (TU] and [TT] the more 
reasonable Bose-Einstein phonon distribution function was used instead of 
the Bose-cone distribution. But the approximate theories developed in these 
papers only allowed us to study phonon evolution which was near to the 
initial state (see [H]). The question of the mesa shaped angular distribution 
did not even occur in [TU] and [TT], because this phenomenon was only fully 
investigated recently [5|. 

In the present paper we rigorously solve the problem of phonon pulse 
expansion into the "superfluid vacuum" of helium. This gives a physical 
explanation of the mesa shape formation and the nontrivial dependencies of 
the mesa width and height on various parameters. 

In Sec. II we obtain the nonlinear equations, which describes the evolution 
of phonon systems created by thermal pulses in superfluid helium. In Sec. Ill 
we find explicitly the family of exact solutions of these equations. They are 
the simple waves, which correspond to the second sound modes in phonon 
systems. In Sec. IV these exact solutions are used to describe the first 
stage of the expansion of a phonon layer in superfluid helium, when only the 
incident waves exist. Here we find the expansion velocity of a phonon pulse 
into the "phonon vacuum" , i.e. into superfluid helium with zero temperature. 
The second stage of the expansion of the phonon layer, when reflected waves 
are formed, are studied in Sec. V. Here we develop an approximate method 
to describe the reflected wave, which allows us to find the average energy 
density and the width of the reflected wave. In Sec. VI by using the theory 
presented in this paper we give qualitative explanations of some experimental 
data from Ref. ([5]). Conclusions are drawn in Sec. VII. 
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2 Equations for the evolution of a phonon 
system in the hydrodynamic approxima- 
tion 

To describe the evolution of a phonon system in superfluid hehum, we start 
from the well-known two-fluid hydrodynamic equations in the non-dissipative 
approximation [12l [13] : 

dp 

— + div{pn^n + Ps^s) = 0; (1) 
at 

— + d^v{S^r^) = 0; (2) 

--^ + V/i+(v,VK = 0; (3) 
at 

dAi , , dT ^ dvnk /.X 

Here p is the density of helium; p„ and ps are densities of the normal and 
superfluid components respectively; Vn and Vg are velocities of normal and 
superfluid components respectively; S is the entropy of unit of volume; p is 
the chemical potential of unit of mass of helium; A = pnw/S; w = v„ — is 
the relative velocity of the normal and superfluid components; T is temper- 
ature. Summation over twice repeated indices k is implied in Eq. (jlj). Eqs. 
([1]) and ([2]) express the conservation of mass and entropy respectively, Eq. 
([3]) expresses the acceleration of the superfluid and Eq. (jH) is a combination 
of momentum conservation and other conservation laws. 

In the experiments of Refs. phonon systems were created by a 

heater immersed in superfluid helium. The temperature of the liquid helium 
was ~ 50 mK. At this temperature, the ambient thermal excitations can be 
neglected. The power of the applied heat pulses was such that the normal 
density pn was very small, pn <^ p- In this case the superfluid velocity and 
variations of pressure can be neglected when we study phonon propagation 
|14j . Taking into account the approximate relation v„ = w -|- w in Eqs. 
([2]) and (jl]), we obtain the system of equations which describes the evolution 
of a phonon system: 

+ *:i,(Sw) = 0; (5) 
g + (wV)) (^w) ^ -VT - gv.^. (6) 



4 



Relations and are the system of four equations for four variables, 
for example, for temperature T and the components of relative velocity w. 
Eqs. ([5]) and ([6]) must be completed by the equations of state of superfluid 
helium 5" = S(T, w^), pn = Pn(T, w"^), in which we have neglected the pressure 
dependence in accordance with the assumptions made above. 

The system of Eqs. (E])-® are linearized for small deviations of the vari- 
ables. They define the propagation of second sound waves and the transverse 
wave, at arbitrary values of the equilibrium relative velocity w, when p„ <^ p. 
The dispersion relation for this case is studied in [13]. 

For a phonon system with a linear energy-momentum relation e = cp, 
where = {dP) / {dp) and c is the first sound velocity of helium, we have 

m- 

Pn _ T _ 27r^ k%T^ 

S C^-W^' 45 n3c3(l-w2/c2)2- 

The system of equations ([S])-®, together with the relations ([7]), describe 
the evolution of phonon systems propagating in superfluid helium. 

The main feature of phonon systems, created in experiments [I]- [5], is the 
high value of the relative velocity w, which has a value close to the first sound 
velocity c. At the same time the temperature of such phonon systems is very 
small, so the normal density satisfies the strong inequality p„ <^ p. As shown 
in [12], using the general conditions of stability for superfluid helium [8], such 
phonon systems are thermodynamically stable, even for values of the relative 
velocity which approach the first sound velocity c, if the temperature is low 
enough. 

In a recent paper [5] detailed measurements of angular distributions of en- 
ergy in phonon systems were reported for different distances from the heater, 
and different heater sizes and powers. Fig{T] illustrates the main idea of the 
experiments [5]. The phonon system (shaded region) is created by a current 
pulse in the heater H, which is immersed in superfluid helium. 

The created phonon system moves in the direction normal to the heater 
(axis z) with a velocity close to the first sound velocity c. The pulse expands 
transversely. This expansion along with the initial phonon energy distribu- 
tion near the heater, determines the angular distribution of energy on the 
detector. By changing the heater size, one can vary the initial transverse size 
Lq of the phonon system, and by the changing heater pulse duration tp one 
can vary the characteristic length dp of the system. 

We consider the simple case when the dependence on z can be neglected 
in the initial value problem for Eqs. ([5]), ([6]). That is valid for sufficiently 
long phonon pulses. To study the behaviour of the transverse expansion of 
phonon pulses, we consider the 2 dimensional case, where all values depend 
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normal fluid moving with v^^w,. 



w , T 

0' 



ct 



H superfluid helium, v =0, T=0 



Figure 1: A phonon system (shaded region) with characteristic temperature 
To and the relative velocity wq in superfluid helium with = and temper- 
ature T = 0. The heater H has a width 2Lq. The characteristic length dp 
of the system depends on the the pulse duration tp. The coordinate frame 
is defined with the ^-axis parallel to the normal to the heater H. D is the 
detector. 

only on one spatial cartesian coordinate, x. The relative velocity w lies in 
the plane xz, i.e. w = {w^, 0, w^). 

Substituting Eqs. ([7]) into Eqs. ([S]), ([1]), and introducing the dimension- 
less variables 

T 

9 = In- Ux = Wx/c, = {wl + wl)/c^, (8) 

I — 

we can rewrite the equations in the matrix form 

jj+M(«..«^)|(^«° 1=0, (9) 

where the matrix M{ux, v^) is equal to 

(5 — ?>V?)Ux 1 — V? —Ux 

^""^ ^ 6(1 -^2)2^, _2u2(1-m2) 2u'ux 

(10) 

We have introduced the dimensionless relative velocity u = w/c in Eqs. 

(ED-dToD. 
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We see that Eqs. flUl)- ffTUl) is a complicated nonlinear system of first order 
equations in partial derivatives, with coefficients that depend explicitly on 
the variable and linearly and quadratically on the variable Ux- We show 
below that Eqs. (19|)- (|T0|) can be solved exactly for arbitrary values of the 
relative velocity w. 



3 Second sound simple waves in phonon sys- 
tems 

In this section we will obtain the class of exact solutions for the system 
of Eqs. fP)- ffTU]) . In Sec. 4 we apply the solutions, along with the initial 
conditions, to the problem of phonon pulse expansion in superfluid helium 
with zero temperature. 

Let us suppose that all three variables of the system of Eqs. (l9|)- (fT0|) 
depend on one unknown function, for example, u = i>{x,t). In other words 
we seek solution in the form oi Q = Q^u), = u^iy), and = v^iy). 
Such solutions are called simple waves (see, for example, [IS]). Substituting 
these relations into equations ([9]), and dividing the equations ([9j) by dvjdx 
we obtain the following matrix equation 




cdt " 



''X 



^ du 



U,. . 11 



We see from Eq. ( ITT]) that the vector with components {dQ / du, dux/ du, dv^ / du) 
is an eigenvector of the matrix M, which corresponds to the eigenvalue 
— ^/|^. Thus, if we know some eigenvector r = {ri,r2,rs) of the matrix 
M (ITOi) . which corresponds to the eigenvalue V, i.e. if Yl\=i^ikrk = Vri 
for any index i G (1,2,3), then from Eq. fill I) we obtain the running wave 
equation for the function z/(x, t) 

^-^ + Viu)^-^ = 0, (12) 
cot ox 

which has velocity cV. The system of ordinary differential equations for this 
solution is 

du du du 

where A is some multiplier. Below we will show that matrix M (ITOl) has three 
pairwise different eigenvalues. Therefore the eigenvectors, which correspond 
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to the same eigenvalue, are coUinear, and A is a coefficient of proportion- 
ality between a given eigenvector r = (ri,r2,r3) and the eigenvector with 
components {dQ / du, dux/ du, du^/ du) . 

Eliminating A from the system (fT3l) , we get the system of ordinary differ- 
ential equations 

dQ dUx du^ , 

— = — = , 14 

r\ ra rg 

where r = (ri,r2,r3) is a eigenvector of matrix ffTOj) of system ([9]), which 
corresponds to the eigenvalue V . The equations ([T^ determine the functional 
dependence between the variables 0, u^^ and in the corresponding simple 
wave solution of system ([9]), (fTOl) . It follows from Eqs. (|T2|) and (fT3|) . that 
any of the variables 9, Ux^ and satisfy the same running wave equation 
(fT2|) as the function i/(a;, t). 

The three eigenvalues of matrix fllUj] are calculated to be 



= ' 

= , (16) 

% = Ux. (17) 

It should be noted that the eigenvalues (|T5l) - (|T7l) of matrix (ITOj) of system 
of Eqs. Q are real, if m < 1, and Vi < V3 < V2. Thus the system of 
Eqs. ([9]) are hyperbolic. In this connection it is interesting to note, that 
the condition w < c, which guaranties hyperbolicity of the system of Eqs. 
©, din]) coincides here with the condition of thermodynamic stability for a 
phonon system [T5] . 

The expressions f lT5|) - f|T7|) give the simple wave velocities Vi,2,3 = cVi,2,3 
in a phonon system. Two of these velocities coincide, as they should, with 
those of second sound propagation in a phonon system [T3| and the third 
with the velocity of transverse waves [T7j. In this case, the values Ux and 
in Eqs. fll5l)- (|T7j) are the mean constant equilibrium values, from which there 
are small perturbations of temperature and relative velocity, which propagate 
with velocities (|T5l) - (|T7l) . Velocities f|T5|) and (fT6|) correspond to the second 
sound propagation velocities, which were found in [H] for arbitrary values of 
relative velocity w. Particularly, at w = 0, we obtain V^i 2 = =Fc/ -\/3 which is 
the well-known phonon second sound velocity, when p„ ^ p. Specific features 
of second sound in anisotropic phonon systems with 7^ 0, i.e. when the 
system is characterised by a certain direction of the relative velocity w, were 
studied in 
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The simple wave velocity f[T7|) corresponds to the one for the so-called 
transverse wave with the dispersion law uj = kvn ^ kw = kw^- The prop- 
erties of the transverse wave and possibility of realising it in phonon pulses, 
were discussed in |17j . 

Let us consider the simple second sound wave with velocity (ITSl) . The 
corresponding eigenvector of matrix (ITIl is equal to 

2u^ - {l-u'^)R 
'l-u^){u^ + R)R 
2(1 -u'^){?>u^ + u^R)^ 

where we denote 









r = 1 
















R = 



- '"'^ (19) 



Substituting the second and third components of the eigenvector (|T8|) into 
the system of equations ([HD, we obtain the equation 

(20) 



R{u^ + R) 2(3m^ + u^R) ' 

which determines the functional connection between and in the simple 
second sound wave with the velocity V\ f lTSl) . 

Eqs. ( l20|l can be integrated by changing variable to variable R f|T9l) . 
After this substitution we get the differential equation with separated vari- 
ables 

dur 2R^dR 



l-ul (i?2_3)(i?2_i)' 
which can be immediately integrated. Its solution is 



(21) 



where Ci is an integration constant. 

Solving Eq. (fT9!) with respect to m^, we find 

«^(i?) = l-2^il^. (23) 

Relations 0221) and 0231) give, in parametric form, the functional connection 
between the x-component u^-, of the dimensionless relative velocity u, and 
its square m^, which corresponds to velocity Vi, in the simple second sound 
wave fflSll. 
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Substituting the first and second component of tlie eigenvector ([TSjl into 
tlie system of equations ([SD, we obtain tlie otlier equation 



,p _ (2^x. - (1 - u^)R)du^ . . 

{l-u^)iu,, + R)R • ^ ^ 

Expressing u"^ in terms of R, and by using Eq. (|23|) in the factor 
outside dUx in Eq. fl2^ . we get the following equation instead of Eq. 



dQ - = - 25 

Now, using Eq. (12T|) in the right-hand side of Eq. fl25l) . we again get a 
differential equation with separated variables. By integrating this equation, 
and then, taking into account the definition ([H]) of variable B, we obtain an 
expression for the temperature T in the simple second sound wave 




T{R) = C2J (26) 



where C2 is an integration constant. 

Thus the relations (1221) , (1231) , and (I26p express in parametric form the re- 
lationships between the x-component of the dimensionless relative velocity 
u, the square of the dimensionless relative velocity u^, and temperature T, 
in the simple second sound wave, which propagates in "superfluid vacuum" 
of ^He with the velocity Vi f|T5l) . The expression for the velocity Vi f|T5|) can 
be transformed into a more simple form using 



HR) = (27) 
R — Ux 

It follows from Eqs. (ETj) . (12^ . and (12^ . that and increase monoton- 
ically when R increases, and T decreases monotonically when R increases. 
Thus starting from any of the variables Ux, v?-, or T at the initial moment of 
time we can determine the corresponding initial function R{x,t = 0). The 
time evolution of the value of R is determined by the running wave equation: 

^ + HW^^O. (28) 

where Vi{R) is the velocity (127|) . The solution of Eq. (l28i) is well-known 

x-cVi{R)t = f{R), (29) 
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where f{R) is an arbitrary function. 

The solution fl2U]) shows that every value of R runs with its own velocity 
Vi{R) in the simple second sound wave. The relations (129|) . (1271) along with 
the expressions (122|) . (123|) . and (l26l) . express in parametric form the spatial 
dependence of the x-component Wx = cux of the relative velocity w, the 
square of the relative velocity w"^ = c^u"^, and temperature T in the simple 
second sound wave at any moment of time. 

It follows from Eq. fl271) . that dVi{R) / dR > 0. So it is clear that a simple 
second sound wave can, in general, give rise to a break in the continuous 
solution . To determine the location and the velocity of the break we should 
return to the energy and momentum conservation laws, however we will not 
be concerned with that problem here. 

The other simple second sound wave, which has the velocity V2 (see Eq. 
(|T6|) ). can be obtained from the solution (l22l) . (l23l) . and fl26l) . which corre- 
sponds to the velocity Vi. For this purpose let us note that if T{x, t), Wx{x, t) 
and w'^{x,t) is a solution of the system (!8|)- (fT0l) . then T{—x,t), —Wx{—x,t) 
and w'^{—x, t) is also a solution. This transformation maps one simple second 
sound wave to the other. 

4 Expansion of a phonon pulse into the "phonon 
vacuum" of superfluid helium 

To study the main features of the transverse expansion of phonon pulses in 
superfluid helium we will solve the following problem. Let superfluid helium 
with Vs = fill up all space. Let us consider the initial conditions 



for the system of Eqs. (15]). pi and ([TUD. 

At the initial time t = 0, in the superfluid helium at T = 0, there is a 
layer of phonons of width 2Lo, temperature T = Tq, and relative velocity 
w = (0,0, Wq) directed along z-axis (see Figl2]). The phonons in this initial 
layer is the normal fluid, and as the phonon density is very low the normal 
fluid density is very small so pn <C p. Therefore we can neglect the superfluid 
velocity and pressure changes when we study the development of the phonon 
system. Thus the set of Eqs. Q and (|T0|) apply to this phonon system. We 
are interested in finding the development in time of the temperature T and 
the relative velocity w of the phonon system. 




(30) 
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Figure 2: A layer of phonons with temperature Tq and the relative velocity 
wq = WqBz in superfluid helium with = and temperature T = 0. The 
layer width is 2Lq. The x-axis is directed perpendicularly to the plane of the 
layer. 



The problem, formulated above, can be solved partly by using the simple 
second sound wave solution found in the previous section. First of all we 
note that the initial conditions fl30|) satisfy the symmetry transformations 
discussed at the end of section 3. Therefore at any time, the temperature T 
and w"^ are even functions of the x-coordinate, and Wx is an odd function. 
Below we will only discuss the solution for domain x > 0. 

The strong discontinuity, which exists a.t x = Lq in the initial conditions, 
disappears immediately after t = 0, as it is clear from physical considerations. 
The phonon gas expands into the superfluid helium with zero temperature, 
forming the forward front of the outgoing second sound wave (see Figl3]). At 
the same time, the initial perturbation, in the form of a weak discontinuity, 
moves in the region x < Lq forming the rear front of the ingoing second 
sound wave (see FigJ3]). The weak discontinuity is at the point where the 
temperature starts to change from its initial value fl30|l . and it moves to- 
wards the coordinate origin a; = with the velocity cVi{ux = 0, Uq = Wq/c^) 
determined by Eq. f|T5l) . Thus until the time to when the intersection point 
reaches x = 0, there is no length scale in the problem, and the solution must 
be self-similar. The value of to is given by 



Vi{ux = 0,ul = wl/c^)\ "VVT 



2 ' 



(31) 



12 




Figure 3: The dependence of the temperature T on the x-coordinate for the 
initial value Wq = 0.95c at time t = (dashed line), t = to/ 2 (dotted line), 
and t = to (solid line). The arrows point to the location of the fronts of the 
ingoing and outgoing waves at times t = to/2 and t = to- 

A self-similar solution is a particular case of a simple wave. Therefore from 
Eq. ( 1291) . taking into account that at the initial time t = the wavefront is 
at x = Lq, we get f{R) = Lq. Thus the function f{R) is reduced to constant 
for the self-similar solution and 

x = Lo + Vi{R)t (32) 

where Vi{R) = cVi{R) is determined by the relation (1271) . 

For any moment of time < t < to, the wave occupies the domain x G 
{xin{t),Xout(t)), where value Xmit) = Lo + Vi{Ro)t determines the position 
of the front of the ingoing wave, and Xoutit) = Lo + Vi{Rout)t determines the 
position of the front of the outgoing wave, which follows from Eq. fl32p (see 
FigJS]). We denote the yet-unknown limits of the variable R as Ro and Rout, 
i.e. R G {Ro,Rout)- Below we show that Rout = +oo, and that the value Ro 
is determined by Eq. fl35|) coinciding with the value R{ux = 0,Mq = Wq/c^) 
from Eq. ffT9|) . 

To determine the integration constants, which are contained in Eqs. f l2^ 
and fl2Bl) . we should join continuously the simple wave solution (1221) . (1251) and 
(|26|) with the to the initial values (l30l) at R = Rq at the front of the ingoing 
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wave. Thus we get the solution 



Ux{R) 




and the value Rq 



R, = J'j^. (35) 

On the front of the outgoing wave, in accordance with the initial condi- 
tions (130|) the phonon temperature is equal to zero. It follows from Eqs. flM|) 
and (133|) that this can only be satisfied by Rout = +00. 

Thus for any t G (0,to), Eqs. ([33]), ([23]), ([32]), and ([27]) determine 

the desired self-similar solution of the initial value problem fllOl) . and 
(13U]) . and give the values of T, w^, = cWa,, and = c^u^ at any point 

X e {Xin{t),Xout{t)). 

From Eq. fl33]) we find that at the front of the outgoing wave, which 
borders with superfiuid helium at zero temperature, the x-component of the 
dimensionless relative velocity u^^ is equal to 

Kut — 1 ^ _/-Ro~i\/-Ro + v3\ 



/J-, , N Ao„t - 1 f Ro - l\ I Rq + V5\ , , 



From Eq. (1271) it follows that this velocity w^out coincides with the velocity 
of the front of the outgoing wave 

Vi{Rout = +00) = w^out = cutout- (37) 

Also from Eq. (127]) in accordance with Eq. (13T]) for the velocity of the 
front of the ingoing wave, we get 



which coincides with the velocity cVi{ux = 0, Mq = Wq/c^) determined by Eq. 

It is interesting that the relative velocity at the front of the outgoing wave, 
which follows from Eq. ([23]) . is equal to the phonon velocity i.e. w'^{Rout = 



14 




Figure 4: The dependence of x-component of the relative velocity on the 
x-coordinate, for the initial value Wq = 0.95c for times t = to/2 (dotted line), 
and t = tQ (solid line). 

+00) = (?. Whereas the temperature of such phonons is equal to zero. 
Therefore such phonon system remains thermodynamically stable [15]. The 
phonon energy density at the front of the outgoing wave tends to zero. 

Figsini 111 El show the spatial dependence of the temperature T (Fig El ), 
the x-component (FigS]) of the relative velocity w, and the square of the 
relative velocity w"^ (Fig|5j) for times t = (dashed lines on Figs|3]and[5]), t = 
to/2 (dotted lines in Figsl3],|4l and [5]), and for t = to (solid lines in FigsJSHH 
and El . These graphs are calculated from Eqs. (|33|), (151, (|23|), (|32l), and (127]) 
taking into account fl35p for the initial value Wq = 0.95c, which corresponds 
to a strongly anisotropic phonon system. We see that temperature (Fig|2D of 
the phonon pulse, expanding in "phonon vacuum", decreases monotonically, 
but the x-component Wx (FigJl]) of the relative velocity w, and the square of 
the relative velocity w"^ (Fig|5]) increase when x-coordinate increases. 

Let us find all the variables at the point x = Lo- If we substitute x = Lo 
in Eq. ( 1521) . then we get the equation Vi{R*) = 0, where we denote as R* 
the value of R, which corresponds to x = Lq at any time t > 0. Taking into 
account Eq. (127|) . we obtain the equivalent equation R*Ux{R*) = 1. Solving 
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Figure 5: The dependence of the absolute value of the relative velocity w 
on the x-coordinate for the initial value Wq = 0.95c for times t = (dashed 
line), t = to/2 (dotted line), and t = to (solid line). 

this equation taking into account the expression fl33l) . we get 

fl. = V3i±I,-,= (^V^)*. (39) 
1-7 \Ro + V3j \R<,-lJ 

where Rq is determined by Eq. ([35]). 

Substituting the value R*, determined by (l39l) . into Eqs. (133|) . (IMIl . 
and (123|) . we obtain the temperature T{R*), the x-component Wx{R*) of the 
relative velocity w, and the square of the relative velocity w'^{R*) at x = Lq, 
i.e. these are the values at the initial pulse boundary, for any time while the 
solution is self-similar. In particular we get 

MR1 = ^^, ^"^R*) = ^-^2- (40) 

If -R > -R*, then Vi{R) > 0, therefore the wave moves along the positive 
direction of the x axis, and if R < R* then vice versa. 

Fig E] shows the dependence of the expansion velocity w^out, of the phonon 
pulse, on the initial relative velocity Wq, calculated from Eqs. (15^ and (15^ . 
We see that for strongly anisotropic phonon systems, when wq ~ c, the 
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Figure 6: The dependence of the expansion velocity Wrcout of the phonon pulse 
on the initial relative velocity Wq. 

expansion velocity Wxout can be very small, compared to the phonon velocity. 
When Wq -C c, then Wxout ~ c, and phonon system expands nearly with the 
phonon velocity. 

The self-similar solution Eqs. ([3S]), dSl, (ES]), (132]), and ([27D, together 
with the unaffected region which remains constant , determine the desired 
solution of the initial value problem for Eqs. ([9]), ffTOj) . and fl30l) everywhere 
only till moment of time to- At this moment of time the front of the ingoing 
wave, which propagates in the region a; > 0, reaches the coordinate origin 
X = 0, and meets the front of the ingoing wave which propagates in the region 
X < in a symmetric way to the a; > wave. At t > to these two waves 
overlap at the centre region of the phonon pulse. Thus a reflected wave arises 
in the region < x < Xref(t), when t > t^. In the region Xrefif) < x < Xoutif) 
the solution remains self-similar, and is described by the same formulae (l33l) . 
([31]), (125]) . (1521) . and (1271) . The point Xrefif) is determined from the condition 
of a continuous join between the reflected wave and the self-similar solution. 

Because we cannot find the analytical solution for the reflected wave, we 
consider in the next section the approximate method for its description. 
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5 Approximate description of the reflected 
wave 

At the time t = to, when the front of the ingoing wave reaches the coordinate 
origin x = 0, the x-coordinate derivative of the energy density has a maxi- 
mum at X = 0, and the initial plateau on the energy density curve has shrunk 
to zero. After the time t = to, when the reflected wave appears, it transfers 
the energy from the coordinate origin, where the energy density is maximal, 
to the periphery, where the energy density tends to zero. As a result, in the 
region < x < Xrefit) the energy density curve has an approximate plateau, 
where the x-coordinate derivative of the energy density is small compared 
to the one in the region where the reflected wave has not reached. This de- 
pendence of the energy density on the x-coordinate (or on the angle between 
the normal to the heater and the detector) has been experimentally studied 
recently in Ref. [5]. At the present time we cannot find the analytical solu- 
tion for the reflected wave because there are no general methods for solving 
nonlinear equations of the form ([9]) and ffTOjl . However preliminary numerical 
calculations confirm the picture presented here. 

Starting from the condition that the reflected wave and the incident wave 
join continuously, and using the energy conservation law, we can find approx- 
imately the average energy density and the average width of the reflected 
wave. For this purpose let us formulate the energy conservation law for one 
dimensional phonon propagation in superfluid helium. 

For a phonon system with a linear energy-momentum relation e = cp, it 
is easy to obtain the energy density 

30 n?c^{i-wyc^y ' ^ ' 

The energy conservation law can be written as 

where is the phonon energy density flux in the x-direction 

~ 45 h^c%l-wyc^)^' ^ ' 

The energy conservation law (1421) along with the expression for the energy 
density flux fH5]) can be derived directly from the system of equations 
and ffTOj) . which describe the propagation of the phonon system in superfluid 
helium. 
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Figure 7: The solid line ACF presents the spatial dependence of the phonon 
energy density E at some time < t < to- 

When the phonon expansion is described by ([9]) and (|T0|) . together with 
the initial conditions (l30i) . the total energy, which is localized at the half- 
plane X > 0, does not depend on time due to the energy conservation law 
fl42p and the symmetry of the initial conditions fl30p . This energy is equal to 
EqLq, where Eq = E(To,wo) = E{Ro) is the initial phonon energy density 

m- 

On FigIT] the solid line represents the spatial dependence of the phonon 
energy density fHTl) at some time < t < Iq, while the phonon pulse expan- 
sion is described by the self-similar solution fl55]) . (^^, fl2^ . ([S2D, and fl27|) 
(line ACF). The energy conservation law leads to the equality of areas ABC 
and CDF in FigH 

At the time t > to, in the central region < x < Xrefit), there is a 
reflected wave. In the region x > Xref{t) the solution is the self-similar wave. 
To find approximately the average energy density and the average width of 
the reflected wave we formally continue this self-similar solution into the 
region x < Xrefif) until some point Xmif), which is characterised by the 
energy density Em{t). We show this in FiglHl where the self-similar wave is 
drawn at some time t > Iq. This value Em{t) along with Xm{t) is determined 
from the condition that the shaded area on Fig|H]is equal to the initial energy 
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Figure 8: The solid line presents the spatial dependence of the phonon energy 
density E at some time t > to. The values Xm and shows the approximate 
position and the average energy density in the reflected wave respectively. 

EqLq. Thus we can write the equation for Em and Xm 

EmXm+ E{R)-^dR = EoLo, (44) 

where Xm = x{Rm) (see Eq. (l32l) ) and Em = E{Rm)- 

Substituting in Eq. (jH]) the expression (!32|) for x-coordinate, we get after 
integration by parts 



{E{Rm) — Eq)Lq 
InrV^iRf-W^dR 

The expression 0451) determines approximately the time tm, when the re- 
flected wave reaches the value Rm > Ro- Note that the integral in the 
denominator of Eq. (145|) tends to zero at Rm = Ro- This is a mathemati- 
cal expression of the equality of areas ABC and CDF in FiglT] at any time 
t < to. 

It should be noted that if the energy density gradient was zero in the 
reflected wave region, then we would have the equality Xmit) = Xrefif) and 
Figj8]would show the exact spatial dependence of the energy density, with the 
characteristic plateau in the central region. As it is clear from a physical point 
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Figure 9: The dependence of the approximate width Xm of the reflected wave 
on time t for the initial value of Wq = 0.95c. 

of view and preliminary numerical calculations, in the reflected wave region, 
the energy density gradient is much smaller than the one in the ingoing wave, 
but not equal to zero. Therefore Xm(t) does not coincide with Xref(t). In fact, 
because dE/dx < 0, Xm(t) < Xref(t). But the difference between Xm(t) and 
Xref{t) is small, if the energy density gradient is small in the reflected wave 
region compared to the one in the ingoing wave. So the theory, developed 
above, allows us to find the average energy density and the average width of 
the reflected wave. This width corresponds to the width of the mesa shape 
in the angular distributions observed in [5]. 

Figini shows the dependence of the approximate width Xm of the reflected 
wave on time t, calculated from Eqs. (jlSD, ([33), dH]), (I33D, (131, (E3D and 
(|27|) . In Figj9]we see that this dependence is nearly linear. This corresponds 
to the constant velocity of the front of the outgoing self-similar wave. In the 
experiments [5] it was found that the mesa width increased with the distance 
from the heater which is the same as the dependence of Xm on time in FigJHl 

The average energy density in the reflected wave as a function of 
time t is presented in Fig{TOl It is calculated from Eqs. fH5|) . fHT]) . fl33l) . 
([31]), fl23l) . and fl271) . We see that the approximate average energy density 
in the reflected wave Em decreases monotonically with time because of the 
expansion of the phonon pulse. The dependence of E^ on time (Fig fTUl) is 
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Figure 10: The dependence of the approximate average energy density Em, 
of the reflected wave, on time t for the initial value of Wq = 0.95c. 

similar to that observed in experiments [5] where the mesa height decreased 
with the distance from the heater. 

6 Comparison with the experimental data 

Initially, the phonon pulse created by the heater (see Fig{T]), is comprised 
of low energy phonons (1-phonons), which immediately start to create high 
energy phonons (h-phonons) [i8\ . The h-phonon creation results in spatially 
smoothing the initial hot nonuniform central region of 1-phonon pulse, which 
is injected by heater into superfiuid helium. However the h-phonon creation 
for short pulses occurs mainly near the heater within a distance of several 
millimeters [18]. The h-phonon creation time can be roughly estimated as 
10 /is for the typical experimental heater powers. This time corresponds to 
the phonon travelling a distance 2.4 mm. The initial conditions fl5Ul) for the 
1-phonon system correspond to the end of the h-phonon creation stage of the 
evolution of the phonon system. So the above theory is for the expansion of 
a cool 1-phonon system. 

After the h-phonon creation rate has become negligible, the temperature 
is To and the relative velocity is Wq. The values of Tq and Wq depend on the 
power applied to the heater. Their values for 1-phonon pulses, with a linear 
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Power 


Initial values 


P, mW 


To, K 


wo/c 


3.125 


0.033 


0.955 


6.25 


0.041 


0.947 


12.5 


0.050 


0.938 


25 


0.065 


0.921 



Table 1: The values of temperature To and the relative velocity Wq (divided 
by first sound velocity c = 238 m/s) for the 1-phonon system with different 
values of the heater power P. 

momentum-energy relation, are calculated using the results of [18] for exper- 
imentally applied powers. The results of these calculations are presented in 
Table 1 for the heater powers used in the experiments [5j. The second and 
third columns show the corresponding values of temperature Tq and the rel- 
ative velocity wo for these powers. Using the theory presented in this paper, 
we can find at any time the spatial dependencies of all the parameters of the 
1-phonon system for different powers, which correspond to the values of Tq 
and Wq in Table 1. 

We should note that the phonon pulses created in experiments [5] are not 
long enough to neglect the ^;-coordinate dependence in the phonon evolution 
system ([5]), ([6]). Moreover the experimental conditions are much more close 
to cylindrical symmetry, then to a plane one. Nevertheless we think that 
the qualitative dependencies, obtained here for a plane phonon layer, can be 
compared to the real much more complicated experimental conditions. 

The theory presented here shows that the initial phonon pulse starts to 
expand into the "phonon vacuum". The simple second sound wave appears 
with the outgoing wave front moving into the "phonon vacuum", and the 
front of the ingoing wave moving to the centre at a; = (see FiglS]). During 
this time, the width of initial plateau decreases with time. So at the time to 
(see Eq. (l3Ti) ) there is no mesa. The phonon pulse propagates a distance of 
Iq = cto during this time. For the initial pulse width 2Lo = 1 mm and for 
powers 3.125 mW, 6.25 mW, 12.5 mW, and 25 mW from Ref. 0, we obtain 
the values of Iq of 2.4 mm, 2.3 mm, 2.1 mm, and 1.9 mm respectively. The 
total distance between the heater and the phonon pulse Ih is Ih = lo+h, where 
/c = 2.4 mm is the length after which the h-phonon creation has effectively 
stopped. Note that this schema of taking into account of h-phonon creation 
is rather rough one, especially for high powers. 

In experiments |5j the minimal distance between heater and detector was 
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Figure 11: The dependencies of the full mesa width 2xm in degrees, on 
heater power for the distance 12.3 mm from the heater, for the heater widths 
2Lq = 1 mm (line 1) and for 2Lq = 0.5 mm (line 2). 

8.2 mm. Thus the decrease in the width of the initial plateau is not observed 
in [5]. It would need measurements of the phonon energy angular distri- 
butions at a distance about 3-4 mm from the heater to compare with this 
prediction of the theory that there is no mesa at ~ 4 mm. 

At time t = to the reflected wave appears. From a physical point of view 
it is clear that in the reflected wave, which occurs in the central region of 
the phonon pulse, the energy density gradient is smaller than in the ingoing 
wave. Thus the existence of the reflected wave should lead to a mesa shape 
in the phonon angular distributions at distances larger than ^ 4 mm. 

In Ref.[5| the detailed measurements of 1-phonon angular distributions 
in superfiuid helium showed a distinct mesa shape. The dependence of the 
mesa height and width, on heater dimensions, distance to the detector and 
heater power, were measured [5J. Using the theory presented in this paper 
we now can give qualitative explanations some of these experimental data. 

In FigHT] we show the dependencies of mesa width Xm on power for the 
distance 12.3 mm from the heater for the heater widths 2Lq = 1 mm (line 1) 
and for 2Lo = 0.5 mm (line 2), calculated from the data in Table 1 and Eqs. 
(gS]), ([32]), ([33]), ([27]), ([35]) and([23]). Here and below we take the distances 
and the heater widths which were used in the experiments [5]. We see the 



24 




I I I I I I I I I I I I 

5 10 15 20 25 30 

heater power (mW) 

Figure 12: The dependencies of the mesa height on power for the distance 
12.3 mm from the heater for the heater widths 2Lq = 1 mm (hne 1) and for 
2Lo = 0.5 mm (hne 2). 

monotonic growth of the mesa width Xm, when power increases as it does in 
the experimental data from [5] , although the measured mesa widths are con- 
siderably smaller than those theoretically predicted. The power dependence 
can be explained in the following way. When power increases, the relative 
velocities Wq decrease (Table 1). This results in increasing the absolute value 
of velocity Vi. Thus the time to (see Eq. ( I3T|) ) decreases, and the mesa forms 
earlier. It should be noted that in the experiments [5] for some cases the 
measured mesa widths start to decrease at high power 25 mW (see Fig. 9 
of [3]). Our theory, which describes the expansion of a cool 1-phonon layer, 
cannot explain this effect. 

In Fig{T2] we present the dependencies of the mesa height Em on power 
for the distance 12.3 mm from the heater for the heater widths 2Lq = 1 mm 
(line 1) and for 2Lo = 0.5 mm (line 2), calculated from the data in Table 1 
and Eqs. (USD, HMD, dSlD, (EB, dSS]) andd^ni). The mesa height E^ increases, 
when the power increases, as it does in the measured data [5j. 

While the mesa width Xm do not depend on the initial temperature Tq, 
which follows from Eqs. (H5|) and (132|) . the mesa height Em depends on Tq . 
It turns out that the increase of Tq with power prevails over the decrease the 
relative velocity Wq, and so an increase in power leads to a quicker expansion 
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Figure 13: The angular width of the full mesa 2xm versus the inverse of the 
heater width for the distance 12.3 mm from the heater, for heater powers 
3.125 mW (line 1) and 25 mW (line 2). 

of the phonon pulse and hence to a wider mesa. 

The increase of Xm with time, shown in Figj9] corresponds to an increase 
of the mesa width with distance. This is observed in experiments [5]. The 
decrease of Em with time, shown in FigUHl corresponds to a decrease in the 
mesa height with distance. This too is observed in experiments [5]. 

The theory developed in this paper allows us to qualitatively understand 
why the mesa width increases with decreasing the heater width [5]. A smaller 
heater width {2Lq) results in an earlier formation of the mesa, and this leads 
directly to an increase in the mesa width. In FigJT3]we show the dependencies 
of mesa width Xm on reciprocal heater width l/(2Lo) for the distance 12.3 
mm from the heater for powers 3.125 mW (line 1) and 25 mW (line 2), 
calculated from the data in Table 1 and Eqs. (05]), ([32]), ([33D, ([22]), §B> 
and fl23l) . But we should note that our theory does not reproduce the linear 
dependence of mesa width on reciprocal heater width, which was observed 
in the experiments ^ (see Fig. 10 in [5]). 
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7 Conclusions 



In this paper, starting from the hydrodynamic equations of superfluid hehum, 
we have obtained the system of equations ([5]) and IQ, which describe the 
evolution of a cool phonon systems, created by thermal pulses in superfluid 
helium for the case p„ ^ p (see FigH]). These equations are simplified to 
the case of one spatial dimension (see Eqs. fl51)- ffTUl) . and FigJ2]). The family 
of exact solutions were found in an explicit analytical form. They are the 
simple waves of second sound for phonon systems. The relations between 
temperature T ( 1261) . the x-component fl22l) of the relative velocity w, and 
the square of the relative velocity (j23D are studied for the simple second 
sound waves in phonon systems. These solutions are used to describe the first 
stage of expansion of a phonon layer in superfluid helium, when only simple 
second sound waves exist. Figsl3|, HI and [5] show the spatial dependence of 
temperature T (FigjS]), the x-component Wx (FigS]) of the relative velocity 
w, and the square of the relative velocity w'^ (FiglS]) for time t = (dashed 
lines on FigsJS] and E]), t = to/ 2 (dotted lines in FigsIS], IH andE]), and for 
t = to (sohd lines in FigsIHUH and [5]). 

We have found the velocity of expansion Wxout of a phonon pulse, propa- 
gating initially in z-direction, into the "phonon vacuum", i.e. into superfluid 
helium at zero temperature. The dependence of the expansion velocity Wxout 
of the phonon pulse, on the initial relative velocity Wq, calculated from Eqs. 
fl36|) and fl35l) are presented in FigEl We see that for strongly anisotropic 
phonon systems, when wq ~ c, the expansion velocity Wxout can be small com- 
pared to the phonon velocity. When wq <C c, then Wxout ~ c, and phonon 
system expands nearly with the phonon velocity. 

In the second stage, after the incident wave reaches the centre of the 
phonon layer, a reflected wave appears. We found the time to Eq. ([21]), when 
this starts. The reflected wave transfers energy from the coordinate origin, 
where the energy density is maximal, to the periphery, where the energy 
density tends to zero. It smooths the dependence of the energy density on the 
x-coordinate, in the reflected wave region. The smallness of the x-coordinate 
derivative of the energy density in the reflected wave region, compared with 
the one in the ingoing wave, results in a mesa shape form, which was observed 
in [5]. 

We developed an approximate theory for the average energy density and 
the average width of the reflected wave (see FiglH]). The calculated depen- 
dencies of the mesa height and mesa width on time (FigsM fTU]) . on the heater 
power (Figsim [12]), and on the heater width (Fig JTS]) show partly the same 
qualitative dependencies as the experimental data in [5j, although our theory 
fails to explain all effects observed. But we think that the main cause of the 
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mesa shape appearance in the experiments [3] are the same as in our simple 
model theory. 
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